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ABSTRACT 



The SED of the lensed z = 6.56 galaxy HCM6A be- 
hind the cluster Abell 370 has been analysed. We 
find clear indications for the presence of dust in this 
galaxy, and we estimate the properties of its stellar 
populations (SFR, age, etc.). From its estimated lu- 
minosity, L (1 — 4) X IO^^Lq, this galaxy ranks as 
a luminous infrared galaxy . This case is then used 
to examine the detectability of high-z galaxies with 
Herschel and ALMA. It is evident that with the use 
of strong gravitational lensing SPIRE/Herschel ob- 
servations could provide very interesting information 
on z > 5 galaxies. Strong synergies between ground- 
based near-IR instruments on 8- 10m class telescopes 
and ELTs, Herschel, the JWST, and ALMA can be 
expected for the exploration of the first galaxies in 
the Universe. 

Key words: Galaxies: starburst - Gravitational lens- 
ing - early Universe - Infrared: galaxies - Submil- 
limeter. 



INTRODUCTION 



Considerable advances have been made in the explo- 
ration of the high-z Universe during the last decade, 
starting with the discovery and detailed studies of 
redshift z ~ 3 galaxies (Lyman break galaxies, 
LBGs), mostly from the pioneering work of Steidel 
and collaborators (cf. Steidel et al. 2003), reaching 
over the z ~ 4-5 galaxies from different deep multi- 
wavelength surveys, up to galaxies at z ~ 6-7 close to 
the end of the reionisation epoch of the Universe (e.g. 
Kodaira et al. 2002, Hu et al. 2002, Cuby et al. 2003, 
Kneib et al. 2004, Stanway et al. 2004, Bouwens et 
al. 2004). To extend the present searches beyond 
z > 6.5 and back to ages where the Universe was 
being re- ionized (cf. Fan et al. 2002), it is manda- 
tory to move into the near-IR bands. This logical 
extension should allow to find and study the galax- 



ies up to z 
telescopes. 



10 with ground based instruments and 



Other approaches could also be used to search for 
galaxies during the re-ionisation epoch. Indeed, 
the recognition of a large negative k-correction of 
dust emission in the sub-mm and mm regime which 
largely overcomes the effect of the inverse square law 
and cosmological surface brightness dimming, and 
the existence of numerous dust rich galaxies up to 
fairly high redshift, indicate that searches with high 
resolution instruments as ALMA could locate high-z 
galaxies in an independent way and probably even 
measure their redshift from fine structure lines (see 
e.g. Blain et al. 2000, Guiderdoni et al. 1999). How- 
ever, to be feasible such searches obviously require 
that the target galaxies are chemically evolved and 
contain sufficient dust. Although quasars up to the 
highest redshift currently known have revealed the 
presence of important quantities of dust (e.g. at 
z = 6.42, Walter et al. 2004), little is known about 
dust in galaxies above z> 5. Actually most observa- 
tions of Lyman break galaxies (LBG) at z > 4 seem 
to show relatively blue colors (and bluer ones than 
LBGs at lower z), which possibly indicates less/little 
reddening in these objects (see e.g. Stanway et al. 
2004, Bouwens et al. 2004, review of Schaerer 2004). 
Finding z 3> 4 galaxies with significant amounts of 
dust, and more generally studying the dust proper- 
ties and content in high-z objects is therefore of great 
interest. 



Here we first report on a recent analysis of a lensed 
z = 6.56 starburst galaxy indicating the presence of 
significant dust extinction (Schaerer & Pello 2004). 
We then use this case to illustrate the detectability of 
such objects with Herschel and ALMA, in particular 
in conjunction with strong gravitational lensing. 
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Abell 370 HCM 6A, z = 6.56 (Hu et al. 2002) Abell 370 HCM 6A, z = 6.56 (Hu et al. 2002) 
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Figure 1. Best fits SEDs to the observations of Abell 370 HCM 6A. The red crosses indicate the corresponding 
model broad band fluxes. Solid lines show the best fit for a template from the BC+CWW group, dotted from 
SB+QSO group, and dashed from the S03+ group Left: Observed spectral range. Right: Predicted SED in 
Spitzer/IRAC domain for best fit models. Dashed lines show the bursts from the BCCWW and S03+ template 
groups. The dotted line is the spectrum of SBS 0335-052 from the SB+QSO group with additional Ay = 1. 
The solid lines show best fits for constant star formation using different extinction/attenuation laws (Calzetti 
starburst law versus SMC law). The solid triangles illustrate the IRAC point-source sensitivity (la) for low 
and medium backgrounds excluding "confusion noise". 



2. STELLAR POPULATIONS AND DUST 
IN A LENSED Z = 6.56 STARBURST 
GALAXY 



The lensed z — 6.56 galaxy HCM6A was found by Hu 
et al. (2002) from a narrow-band survey in the field 
of the lensing cluster Abell 370. Its redshift is estab- 
lished from the broad-band SED including a strong 
spectral break, and from the observed asymmetry of 
the detected emission line identified as Lya. 

We have recently analysed the SED of this object by 
means of quantitative SED fitting techniques using 
a modified version of the Hyperz code of Bolzonella 
et al. (2000). The observed VRIZJHK' data are 
taken from Hu et al. (2002). The gravitational mag- 
nification of the source is /i = 4.5 according to Hu et 
al. The main free parameters of the SED modeling 
are the spectral template, extinction, and the red- 
dening law. Empirical and theoretical templates in- 
cluding in particular starbursts and QSOs (SB-I-QSO 
templates), and predictions from synthesis models 
of Bruzual & Chariot (BC+CWW group) and from 
Schaerer (2003, hereafter S03) are used. 

Overall the SED of HCM 6 A (see Fig. 1) is "red- 
dish" , showing an increase of the flux from Z to H 



and even to K^. From this simple fact it is already 
clear qualitatively that one is driven towards stellar 
populations with a) "advanced" age and little ex- 
tinction or b) constant or young star formation plus 
extinction. However, for HCM6A a) can be excluded 
as no Lya emission would be expected in this case. 

Quantitatively, the best solutions obtained for three 
"spectral template groups" are shown in the left 
panel of Fig. 1. The solutions shown correspond to 
bursts of ages ^ 50-130 Myr and little or no extinc- 
tion. However, as just mentioned, solutions lacking 
young (< 10 Myr) massive stars can be excluded 
since Lya emission is observed. The best fit em- 
pirical SB-I-QSO template shown corresponds to the 
spectrum of the H ii galaxy SBS 0335-052 with an 
additional extinction of Ay = 1. On the basis of the 
present observations a narrow line (type II) AGN 
cannot be ruled out. To reconcile the observed SED 
with Lya, a young population e.g. such as SBS 0335- 
052 or constant SF is required. In any of these cases 
fitting the "reddish" SED requires a non negligible 
amount of reddening. 

To illustrate the typical range of possible results we 
show in Fig. 2 a contour map and the correspond- 

^The significance of a change of the SED slope between JH 
and HK seems weak, and difficult to understand. 
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Chi2 contours (z — 6.56 galaxy: Hu et al. 2002) 
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Figure 2. contour plots showing solutions in ex- 
tinction - age diagrams. The best solutions are in- 
dicate by the black dot. Equidistant levels with a 
spacing of 0.5 are shown. The 2D 68% confidence 
region (corresponding to A^^ = 2.3j is delimited by 
the solid thick black line. The (ID) 68 % confidence 
region for Ay (^X^ — ^.t each given age is de- 
limited by the dashed thick black line. Plot for so- 
lutions using a solar metallicity burst template from 
the SOS-h template group and the Calzetti attenua- 
tion law. The solutions indicate an non-negligible 
extinction, but no constraint on age. Discussion in 
text. 



ing confidence intervals for constant star formation, 
solar metallicity models (S03+ template group) and 
reddened with the Calzetti law. We see that for a 
given age Ay is typically ~ 0.5-1.8 mag at the 68 % 
confidence level. For obvious reasons, no constraint 
can be set on the age since the onset of (constant) 
SF. Hence, from the photometry of HCM 6A and 
from the presence of Lya we are led to conclude that 
this object must suffer from reddening with typical 
values of Ay ~ 1. for a Calzetti attenuation law. A 
somewhat smaller extinction {Ay ~ 0.4) can be ob- 
tained if the steeper SMC extinction law of Prevot 
et al. (1984) is adopted. From the present data it 
is not possible to distinguish the different extinc- 
tion/attenuation laws. 

From the best fit constant SF models we deduce an 
extinction corrected star formation rate of the or- 
der of SFR(UV) ~ 11 - 41 Mq yr-i for a Salpeter 
IMF from 1 to 100 Mq or a factor 2.55 higher for 
the often adopted lower mass cut-off of 0.1 Mq. For 
continuous SF over timescales tsp longer than ^ 10 
Myr, the total (bolometric) luminosity output is typ- 



ically ~ lO^"^ Lq per unit SFR (in Mq yr"!) for a 
Salpeter IMF from 1-100 Mq, quite independently 
of metallicity. The total luminosity associated with 
the observed SF is therefore L ~ (1 - 4) x 10^^ Lq, 
in the range of luminous infrared galaxies (LIRG). 
For ^SF ^ 10 Myr the estimated stellar mass is 
« tsF X SFR ~ (1 - 4) X 10« Mq. Other prop- 
erties such as the "Lya transmission" can also be 
estimated from this approach (see Schaerer & Pello 
2004). 

It is interesting to examine the SEDs predicted by 
the various models at longer wavelengths, including 
the rest-frame optical domain, which is potentially 
observable with the sensitive IRAC camera onboard 
the Spitzer Observatory and other future missions. 
In the right panel of Fig. 1 we plot again the 3 best 
fits. We see that these solutions have fiuxes compara- 
ble to or above the detection hmit of IRAC/Spitzer 
^. On the other hand the strongly reddened con- 
stant SF or young burst solutions do not exhibit a 
Balmer break and are hence expected to show fiuxes 
just below the IRAC sensitivity at 3.6 /xm and signifi- 
cantly lower at longer wavelengths. As Lya emission 
is expected only for the reddened SEDs the latter 
solutions are predicted to apply to HCM 6A. If pos- 
sible despite the presence of other nearby sources, 
IRAC/Spitzer observations of HCM 6 A down to the 
detection limit or observations with other future 
satellites could allow to verify our prediction and 
therefore provide an independent (though indirect) 
confirmation of the presence of dust in this high-z 
galaxy. 



3. Z>6 STARBURSTS: WITH HERSCHEL 
AND ALMA, AND NOW . . . 

Let us now assume that starburst galaxies with dust 
exist at z > 6 and briefiy examine their observabil- 
ity with facilities such as Herschel and ALMA. To 
do so we must assume a typical galaxy spectrum in- 
cluding the dust emission. For simplicity we here 
adopt the SED model by Melchior et al. (2001) based 
on PEGASE.2 stellar modeling, on the Desert et al. 
(1990) dust model, and including also synchrotron 
emission. Their predicted SED for a galaxy with an 
SFR and/or total luminosity quite similar to that 
estimated above for HCM6A is shown in Fig. 3. 

Figure 3 shows the exquisite sensitivity of ALMA 
in the various bands allowing in principle an easy 
detection of such objects up to redshift ~ 10 or even 
higher! 

On the other hand, with the sensitivity of PACS and 
SPIRE blank field observations of such an object are 
limited to smaller redshift {z< 1-4). However, al- 
ready with a source magnification of ~ 3-10 or 

^Seejhttp : //ssc . spitzer . caltech. edu/irac/sens .html 
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Moderate starburst - U= 1 .8 x 1 0" U - SFR = 32 M„ yr" 
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Figure 3. Predicted spectrum for a "moderate" star- 
hurst with SFR=32 Mq yr'^ and L - 1.8 x 10"Lq 
placed at redshift z =0.1, 0.5, 1, 2, 3, 5, 10, 20, 
and 30. The thresholds of the JWST (here NGST), 
PACS and SPIRE onboard Herschel, and ALMA are 
also presented. Figure taken from Melchior et al. 
(2001) with kind permission 



more the "template galaxy" shown in Fig. 3 becomes 
observable with SPIRE at ~ 200-670 ^m. In fact, 
such magnifications (and even higher ones) are not 
exceptional in the central parts of massive lensing 
clusters. E.g. in our near-IR search conducted in 
two ISAAC fields 2.5x2.5 arcmin^) of two lensing 
clusters, a fair number (~ 10-20) of z > 6-7 galaxy 
candidates with n>5 are found (Pello et al. 2004; 
Richard et al. 2004, in preparation). More than half 
of them have actually magnifications > 10! Such 
simple estimates show already quite clearly the po- 
tential of strong gravitational lensing to extend the 
horizon of SPIRE/Herschel observations beyond red- 
shift z>5! 

Obviously a more rigorous feasibility study must also 
address the following issues: How frequent is dust 
present in high-z galaxies? and up to what redshift? 
We now have some indications for dust in one lensed 
z = 6.56 galaxy (see Section 2) and of course in 
high-z quasars. But how general/frequent is this? 
How typical is the SED adopted above? The long 
wavelength emission due to dust depends on vari- 
ous parameters such as metallicity, the dust/gas ra- 
tio, geometry, the ISM pressure etc. Furthermore 
spatial resolution and source confusion are key is- 
sues which must be addressed and which should vary 
quite strongly between blank fields and cluster envi- 
ronments. Last, but not least, the field of view of 



the various instruments is determinant for the effi- 
ciency with which high-z candidates can be found 
and studied. Several of these isssues have already 
been partly addressed earlier (cf. the 2000 Herschel 
conference proceedings of Pilbratt et al. 2001, also 
Blain et al. 2002). 

It is evident that various ground-based and space 
bourne facilities and instruments will be used to- 
gether to provide an optimal coverage in wavelength, 
spatial resolution and field size, and to obtain imag- 
ing as well as spectroscopy. Near-IR wide field im- 
agers and near-IR multi-object spectrographs on 8- 
10m class telescopes and later with ELTs will un- 
doubtably "team up" with the JWST, Herschel and 
ALMA to explore the first galaxies in the Universe 
and their evolution from the Dark Ages to Cosmic 
Reionisation. 
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